Kinetic energy of a trapped Fermi gas interacting with a Bose-Einstein condensate 
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We study a confined mixture of bosons and fermions in the regime of quantal degeneracy, with 
particular attention to the effects of the interactions on the kinetic energy of the fermionic compo- 
nent. We are able to explore a wide region of system parameters by identifying two scaling variables 
which completely determine its state at low temperature. These are the ratio of the boson-fermion 
and boson-boson interaction strengths and the ratio of the radii of the two clouds. We find that the 
effect of the interactions can be sizeable for reasonable choices of the parameters and that its exper- 
imental study can be used to infer the sign of the boson-fermion scattering length. The interplay 
between interactions and thermal effects in the fermionic kinetic energy is also discussed. 
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I. INTRODUCTION 

The achievement of Bose-Einstein condensation in trapped gases has opened new opportunities for investigating 
the low temperature behaviour of dilute quantum systems. Recent experimental studies have been also addressed to 
the search of degeneracy effects in Fermi gases ^|,[5j and in mixtures of bosons and fermions g| . First experimental 
evidences of these effects have been recently reported in j?]] . Due to the Pauli exclusion principle the effects of the 
interactions in a Fermi gas are much weaker than for a Bose condensate |g|-|l3| . Therefore, the kinetic energy dominates 
the behaviour of the Fermi gas and is a clear indicator of its quantum degeneracy. 

Interest in boson-fermion mixtures is stimulated by the fact that, whereas in a one-component spin-polarized Fermi 
gas the absence of interactions leads to long thermalization times which hamper the process of evaporative cooling, 
the collisions between the two species in a mixture can ensure fast thermalization (the so-colled sympatethic cooling 
jl4]-|l6|] ) . The kinetic energy of the Fermi component in such a mixture could be measured by time-of- flight techniques 
in an ideal experiment in which the confining potential is suddenly switched off after a fast expulsion [i.e. on a time 
scale shorter than the boson-fermion collision time) of the bosons from the trap. In this way one can avoid the effects 
of the interactions during the expansion of the gas. Alternatively, the kinetic energy could be obtained from inelastic 
photon scattering at high momentum transfer as recently shown in the case of a trapped Bose gas fl7|| . 

In this paper we analyze the behaviour of the kinetic energy of the fermionic component in a boson-fermion mixture. 
We find that the kinetic energy can be significantly affected by the interactions of the Fermi gas with the Bose-Einstein 
condensed cloud for reasonable choices of the system parameters. An important consequence of this purely quantum 
effect is that one can measure the sign and the strength of the boson-fermion scattering length by using the Bose 
component as a tunable device to change the effective potential felt by the fermions [Q. We assume a positive 
boson-boson scattering length, with a view to applications to mixtures of Rb-K and Na-K. 



II. INTERACTING FERMI-BOSE MIXTURES 



For the description of the mixture at finite temperature we adopt the semiclassical three-fluid model already 
developed in ref. p9[ |. We consider a system of Nf fermions of mass mj and Nb bosons of mass mj, confined by 
external potentials Vj^(r) = ^mf i ) uj'j b r 2 with frequencies Uf and u>b- The external potentials are assumed to be 
spherically symmetric, the asymmetric case requiring simply a change of variables in the framework of the semiclassical 
approximation that we adopt. We include the interaction between bosons through the scattering length abb and the 
interaction between bosons and fermions through the scattering length abf. By assuming that a single spin state is 
trapped for each component of the mixture, we can safely neglect the fermion-fermion interaction which is inhibited 
by the Pauli exclusion principle. Extensions to multi-spin configurations can be naturally made within the present 
formalism. In the following we will neglect the possibility of a superfiuid phase for the fermionic component (for 
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a discussion of the BCS transition in trapped Fermi gases see |2(J,|2lj ) as well as the possibility of the expulsion of 
the bosons from the center of the cloud (this phase separation is expected to occur only for values of the coupling 
strengths such that the mean field contribution of the fermions is comparable with the mean field contribution of the 
bosons. This requires values of the parameters such that rifa 3 b ~ 1. The resulting system is no more a dilute one, see 

©ID- 

The spatial densities of the condensed bosons (n c ), of the bosonic thermal component (n nc ) and of the fermions 
(nf) are determined by the self-consistent solution of the following equations: 

n c (r) = — (nt, - VXt(r) ~ 2g bb n nc {r) - 9bfn f (r)) , (1) 



n, 



9bb 



and 



Here, the effective potentials acting on the thermal boson cloud and on the fermions are given by 

7 eff 

and 



V eff( r ) = V L(r) + 2g bb n c {r) + 2g bb n nc (r) + g bf n f (r) (4) 



= V Lti r ) + 9bfn c {r) + g bf n nc (r), (5) 

where we have introduced the notations (3 = (KbT)' 1 , g bb = Anh 2 a bb / m b and g b f = 2nh 2 a b f /m r with m" 1 = m^ 1 + 
mj 1 . The chemical potentials fi b and /if are determined by the normalization conditions N b = f (n c (r) + n nc {r)) d 3 r 
and Nf — J nAr) cPr, which ensure the self-consistent closure of the model. 

Equations (0-p|) have been derived from a grand-canonical Hamiltonian in which the interactions are included in 
a mean- field Hartree-Fock approximation j24] , by employing the semiclassical approximation for the bosonic thermal 
cloud and for the fermions and by taking the strong coupling limit N b a bb /a,} lo ^> 1 for the wave- function of the 
condensate, with aho = (h/mhuit,) 1 ^ 2 the bosonic harmonic oscillator length. Upon averaging the Hamiltonian on the 
equilibrium state of the system at finite temperature we obtain the energy as the sum of various contributions: the 
kinetic and the external confinement energy for each of the species, as well as the boson-boson and boson-fcrmion 
interaction terms. One has 

E = E Ln + E Lt + E kin + E ext + E int + E \L 



(S) 3/2r5/2 /^^^ + K^) 3/2r5/2 /^ ( - )dV+ 

VL(r) (n c (r) + n nc {r)) d\ + J V e Ur)n f (r) d\ + 

/ ( n cM + 2< c (r)+4n c (7>„ c M) d 3 r + g bf J (n c (r) + n nc (r)) n f (r) d 3 r (6) 

where f p (z) = Tip)- 1 J y^' 1 dy / '(z" 1 e y + 1), g p {z) = Tip)- 1 J yP- 1 dy / '(z^e" ~ 1) are the usual Fermi and Bose 
functions and Zf tb — exp(/3(^/.b — V/ff(r))) . The release energy is obtained by setting the confinement potentials 

in eq. (||) to zero. This quantity can be measured via time-of-fiight experiments. 
At low temperature the thermal component n nc can be safely neglected in the right hand side of equations (0),(^) 
and (0). Similarly, when the fermionic density n/ is much smaller than the density n c of the Bose condensate, its 
contribution in the right hand side of the same equations can be dropped. This is valid, for example, if the interaction 
strengths g b f and g b b have comparable size and if the trapping potential v/ xt is not too stiff with respect to V^ xt . In 
this case the density profile of the condensate is not affected by the interaction with the fermions and the effective 
potential felt by the fermions takes the simplified form: 

. ~m/cj?(l - 7 )r 2 + 9 -^-ii h for r < R b 
Vi ff {r) = {\ 9bb (7) 



?r 2 for r > R, 
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where 



7 =^^4 (8) 
g bb mjujj 

and Rb = {2^/mbUi 2 ) 1 ^ 2 is the radius of the condensate cloud. The potential (0) depends on temperature through the 
boson chemical potential /if,, which determines the radius Rb- Full numerical calculations, including the contribution of 
the thermal Bose and Fermi components, show that this simplified model (hereafter called the double-parabola model) 
describes very well the main features of the system below the critical temperature for Bose-Einstein condensation. 

In the double-parabola model, if the parameter 7 in eq. (|^) is smaller than unity, the potential felt by the fermions 
has its minimum in the center of the trap. In this situation two limiting cases can be envisaged by comparing Rb with 
the radius of the Fermi cloud. This is approximatively given by the Fermi radius Rf = (2Ep /rrifui 2 ) 1 / 2 calculated 
in the absence of interactions, where Ep = (pNfY^ftiJf is the Fermi energy. In the limit Rb <€. Rf the number of 
bosons is much smaller than the number of fermions and thus the interactions play a minor role. Instead, in the limit 
Rb 3> Rf the fermionic cloud feels a harmonic trapping potential with a renormalized frequency uif — u>f(l — r )) 1 ^ 2 - 
Finally in the case 7 > 1 the repulsive interaction with the bosons is stronger than the external potential. The 
effective potential M) then exhibits a local maximum at the center of the trap. 



III. SCALING AND ROLE OF THE INTERACTIONS AT T — 



Let us begin our discussion in the framework of the double-parabola model introduced in the previous section. In 
the case 7 < 1 we can give a simple approximate solution of the model by a variational minimization of the energy 
functional in eq. (0). At T = this reads 



E(T = 0) = jj^-( 67 r 2 ) 2 / 3 J nf{r) d 3 r+ \ \ Z , ,,)„,.„■, ,/ V + 



V e f xt (r)n f (r) d 3 r+'^j n 2 c (r) d 3 r + g bf J n c (r)n f (r) d 3 r. (9) 

The variational approach (see the details in the Appendix) explicitly shows that the relevant properties of the system 
depend on the various parameters of the model through two adimensional combinations, which are the parameter 7 
in eq. (|^) and a parameter x given by 

( I5a bb N b \ 1/5 




2m b uj b V a ho (6N f ) 5 / 6 



At given 7, the ratio of the sizes of the two clouds in the absence of interactions determines the deviation of the 
kinetic energy of the Fermi component from its ideal-gas value. 

We have checked numerically that the scaling in these two variables is satisfied with good accuracy also by the 
full numerical solution of eqs. (|l]-||) at zero temperature for any value of 7. The description of the system with only 
two scaling parameters instead of the eight original ones entering eqs. (|]-||) represents a major simplification of the 
problem. In view of this property, in the following we shall present our discussion in terms of the scaling parameters 
x and 7. 

In fig. |l].a we show a plot of the kinetic energy as a function of x at zero temperature for different values of 
7 < 1. As x increases, the kinetic energy of the fermions goes from its non-interacting value E° in = 3NfEp/8 to the 
strong-coupling limit E^m = ^k„(1 — t) 1 ^ 2 - As a first result of our analysis we see from fig. [l].a that there is a clear 
correspondence, for a fixed value of x, between the value of the kinetic energy and the value of 7. Therefore the sign 
and the strength of the ratio between the boson-fermion and boson-boson coupling constants could be inferred from 
a measurement of the fermion kinetic energy. 

In the case 7 > 1 the fermions are expelled from the center of the trap |l^,^2| and form a shell around the bosons as 
Nb increases with respect to Nf. In this case the kinetic energy of the fermions (fig. [l].b) tends to zero when x — > 00. 

For completeness we have also analyzed the behaviour of the mean square radius of the fermionic cloud as a function 
of x. For 7 < 1 the asymptotic value at large x is larger (smaller) than in the ideal case for 7 > 0(< 0) (see fig. |].a). 
For 7 > 1 the mean square radius increases indefinitely with increasing x (fig. |^.b). These behaviours are immediately 
understood in terms of the behaviour of the kinetic energy shown in fig. \U. 
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IV. THE ROLE OF TEMPERATURE 



Let us finally examine the temperature dependence of the kinetic energy of the Fermi gas. To this purpose we have 
solved self-consistently the full set of eqs. Of course, in the classical regime the kinetic energy is insensitive to 

interactions. As quantum degeneracy sets in at T < Tp, where Tp = Ep/Kb is the Fermi temperature, deviations 
from the classical value 7>N $KbT j2 become apparent. We show in fig. || the predicted behaviour for a given choice of 
the parameters of the mixture. The role of the interactions decreases as temperature increases. This can also be seen 
in fig. U where we plot the kinetic energy as a function of x for a choice of different temperatures. 

The scaling behaviour described in sec. |l| is less accurate at finite temperature. This is easily understood from 
the fact that the approximations leading to eq. (pi) become less justified as temperature increases. 



V. CONCLUSIONS 



We have presented a broad study of the kinetic energy of the fermionic component in a mixture of bosons and 
fermions in the so-called Thomas-Fermi regime {Nbdbb/O'ho 3> 1). We have shown that at zero temperature the kinetic 
energy, as well as the mean square radius of the Fermi component, exhibit an important scaling behaviour in the 
relevant parameters 7 and x defined in eqs. (||,[l0|). This has allowed us to give a systematic investigation of these 
physical properties on a wide range of system parameters and to understand the role of the interactions between the 
Fermi gas and the Bose condensate. 

In particular, we have found that the shift of the fermionic kinetic energy due to the interactions becomes sizeable 
at appreciable values of the parameter x mesauring the relative radii of the two clouds. This effect could be used to 
infer the sign of the boson-fermion scattering length from measurements of the fermion kinetic energy. 

Finally, the role of temperature has been investigated within the self-consistent numerical solution of the full set 
of equations for the coupled boson-fermion mixture and the interplay between thermal and interaction effects on the 
kinetic energy has been demonstrated. 

This work is supported by the Istituto Nazionale per la Fisica della Materia through the Advanced Research Project 
on BEC. One of us (L.V.) acknowledges the hospitality of the Scuola Normale Superiore di Pisa during part of this 
work. 



APPENDIX A: VARIATIONAL MODEL 



The scaling behaviour discussed in section [III| can be explicitly predicted by a variational approach which turns out 
to be very accurate in reproducing the numerical results at T — in the case 7 < 1. The variational method is based 
on the minimization of the energy functional (|^) within a restricted class of functions. 

For 7 < 1 it is convenient to describe the fermionic cloud as if it were embedded in an effective potential V var (r) = 
|mjwj ar r 2 , where the frequency oj var is a variational parameter. The corresponding fermionic density profile is 



n/(r) = l/{6n 2 )(2m f /h 2 ) 3 / 2 (E 
the form 



var 
F 



l f UJ va 



r- 2 /2) 3/2 and the expression for the variational energy functional takes 



E, 



h o 



-NfEj 



ul,. 



fjbf 
9bb 



8 " uj f 

min(i?p ar ,i? b ) ^ 

6^2 



N f Ep 



Uf 



2m f 



3/2 



3/2 



Hb - -m b u? h r 2 ) d 3 r. (Al) 



Here E F ar = (6iV/) 1 / 3 ?ia; t , ar and R F ar = (2E F ar /m/w^) 1 / 2 are the the Fermi energy and the Fermi radius calculated 
with the frequency uj var . The bosons are described by the Thomas-Fermi inverted parabola n b {r) = g^, — 
m b u> 2 r 2 /2) with [i h = huj b (15N b a b b/a ho ) 2/5 /2. 

The integral in eq. (Al)can be carried out analytically, with the result 



E var (x,j, a) = -NfEp x < 



x 2 


x 2 


< 




x z 


X 


a 2 " 





for a < 1 
for a > 1 



(A2) 
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where a 2 = x 2 uj var /ujf and 



P(a) 



2 r 

9tt 



tyjl - a 2 {9 - 18a 2 



40a 4 



16a 6 ) + 3(-3 + 8a 2 ) arcsin(a) 



(A3) 



HW( 15 «* Nb Y /5 . (A4) 



2m b uj b \ a ho (67V/) 5 / 6 J 

The condition dE var /da — determines the value of a and hence of w var . This equation has to be solved 
numerically, except for a > 1 where the model gives the result u> var = Cijt. 



The expression (A2) allows an explicit identification of the scaling variables introduced in sec. III. In fact, the 
quantity E var /NfEp at its minimum depends only on x and 7. 

Of course the variational estimate gives an upper bound for the total energy. This bound is very close to the value 
obtained by solving the Schrodinger equation with the potential (^). Typical deviations are less than 1 % of the 
energy. 
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FIG. 1. Kinetic energy of the Fermi component in units of E® in — 3NfEp/8 as a function of the scaling parameter x for 
different values of 7. a): for 7 < 1 taking the following values: 7 = 0.9 (short dahses), 7 = 0.5 (dots), 7 = 0.1 (solid), 7 = —0.1 
(long dashes), 7 = —0.5 (dot-short dashes) and 7 = —0.9 (dot-long dashes), b): for 7 > 1 ranging from 7 = 1.1 (solid curve) 
to 7 = 1.7 (long dashed curve) in steps of 0.2. 
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FIG. 2. Mean square radius of the Fermi component in units of jj 2 )o = 3NfR%/8 as a function of the scaling parameter x 
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for different values of 7. a): same values and notations as in figure hla except for the value 7 = 0.7 (long dashes-short dashes) 
which replaces the value 7 = 0.9 of figure 0a. b): for 7 > 1, ranging from 7 = 1.1 (solid curve) to 7 = 1.7 (long dashed curve) 
in steps of 0.2. 
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FIG. 3. Kinetic energy of the Fermi component as a function of the reduced temperature T/Tf for a choice of the scaling 
parameters corresponding to m/ = mb = 39 a.u., u>f = uib = 2n x 100 sec -1 , Nf = 10 4 , Nb = 10 6 , abb = 92 ao and at/ = 46 ao 
where ao is the Bohr radius. The solid curve refers to the interacting system and the dashed line to the non-interacting one. 
The dotted curve shows the classical result. 
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FIG. 4. Kinetic energy of the Fermi component as a function of the scaling parameter x for 7 = 0.5 at different temperatures, 
as indicated in the figure. 
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